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ABSTRACT 

We report on the Swift discovery of a second high-amplitude (factor 100) outburst of the Seyfert 
1.9 galaxy IC 3599, and discuss implications for outburst scenarios. Swift detected this active galactic 
nucleus (AGN) again in February 2010 in X-rays at a level of (1.50±0.11)x 10^® W (0.2-2.0 keV), 
which is nearly as luminous as the first outburst detected with ROSAT in 1990. Optical data from 
the Catalina sky survey show that the optical emission was already bright two years before the Swift 
X-ray high-state. Our new Swift observations performed between 2013 and 2015 show that IC 3599 is 
currently again in a very low X-ray flux state. This repeat optical and X-ray outburst, and the long 
optical duration, suggest that IC 3599 is likely not a tidal disruption event (TDE). Instead, variants of 
AGN-related variability are explored. The data are consistent with an accretion disk instability around 
a black hole of mass on the order 10®-10^ Mq; a value estimated using several different methods. 

Subject headings: galaxies: active — galaxies: individual (IC 3599) — galaxies: Seyfert 


1. INTRODUCTION 

The highest amplitudes of X-ray variability recorded so 
far among galaxies have been detected from the cores of a 
few quiescent galaxies, and have been interpreted as the 
tidal disruption and accretion of stars b y dormant super - 
massive black holes (SMBHs; review bv lKomossall2012l) . 
These events happen when a star gets too close to the 
central SMBH and is disrupted due to tidal forces. Part 
of the stellar debris will form an accretion d isk and will 
be swallowed by the SMBH le.g.. lRee'slll988D . producing 
a luminous accretion flare. One of the most spectacu¬ 
lar of these events, S wift J1644-I-57' was discovered in 
Marc h 2011 by Swift ([Bloo^^rami^mjBunowsi^t^I] 
1201 It) . Before Swift (jCehrels et a l. I |2004[) . one of 


the most successful surveys to discover TDEs was the 
ROSAT All-Sky S urvey (RASS), during which several 


events were found I 

Bade et al.l 119961: IGruoe et al. 

|1999|; 

IKomossa & Greiner 

119991: Komossa & Bade 111990 

). 


The Seyfert 1.9 galaxy IC 3599 (Zw 159.034; IRXS 
J123741.2-b264227; azooo = 12’"37'"41."2, ^2000 = 
-b26°42'27". z=0.0215: ICnipe et ahlllMa lBrandt7et all 
119951 : iKomossa fc Bade 111999( 1 was noticed in the RASS 
due to its very bright X-ray emission. However, when 
ROSAT re-observed IC 3599 again about a year later its 
flu x (0.1-2.4 keV) had dropped by a fact or of more than 
60 (|GruDe et al.ll99dlBrandt et al.lll995D . and continued 
to fade, implying a total amplitude of variability larger 
than a factor 100. Its X-ray spectrum was very soft, 
and the high-state X-ray emission was accompanied by 
stron g, broad emission lin es from hydrogen, helium, and 
iron (iBrandt et al.l[T99^. which faded significantly in 
subse quent years (|GruDe et al.lllOO'st IKomossa fc Badel 
II 999 II . This kind of emission-line flaring is very rare, and 
has only recently been observed in a few other galax- 
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ies fe.g.. IKomossa et al. 11200^ . A broad Ha line com¬ 
ponent was still present 7 years after the high-state of 
IC3599 (and still is in a recent SDSS spectrum; see 
below), which led to a Seyfert 1.9 classification of the 
galaxy (IKomossa fc Bade~l 1199911 . The classification of 
IC 3599 as an AGN is consistent with its radio detection 
(iBower et al.l I2013I1 a nd with its AGN-dominated MIR 
spectrum (iSani et al.ll2010H . The cause for the outburst 
of IC 3599 remained unknown, and several possibilities 
were considered: high-amplitude Narrow-line Seyfert 1 
(NLSl) variabilitjQ, an accretion disk instability, or a 
TDE. Its bright high-ionization iron lines were success¬ 
fully modeled by photoionization of high-density gas with 
the properties of a CLR (coronal line region), illumi- 
nated by a flare with a s trong EUV - soft X-ray bump 
(|Komossa fc BadFI 11999( 1 . The faint coronal lines still 
present years after the outburst did not vary further, in¬ 
dicating that they are permanent, and similar to those 
seen in other Seyfert galaxies. 

When Swift observed IC 3599 for the first time in 2010 
February, it detected a new outburst. The luminosity 
level of 1.5x10®® W (0.2-2.0 keV band) was only a factor 
of 2 below t he value seen during the RASS (|Grupe et al.l 
I1995L1200111 . After we noticed the high-state, we trig¬ 
gered new Swift follow-ups of IC 3599 which are reported 
here, along with other multi-wavelength data and a re¬ 
discussion of outburst scenarios. Results of our study 
were first pre sented bv IKomossa et ^ (12015( 1. and inde¬ 
pendently bv ICampana et ahT i 2015(1 . who explored the 
possibility that IC 3599 is experiencing repeat tidal strip¬ 
ping of an orbiting star. 

Throughout this paper spectral indices are denoted 
as Eh oc v~°‘. Luminosities are calculated assuming 
a ACDM cosmology with r2M=0.286, Ha=0.714 and a 
Hubble constant of Hq=70 km s“® Mpc“^. This results 
in a luminosity d istances £>=91 .4 Mpc using the cosmol¬ 
ogy calculator bv I Wright ((2006(1 . All errors are la unless 


^ Note, however, that IC 3599 is very likely not a NLSl galaxy, 
since it does not show any Fell emission, neither at outburst nor 
quiescence, and since its Seyfert 1.9 classification indicates that the 
inner, high-velocity BLR is obscured 
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stated otherwise. 

2. OBSERVATIONS AND DATA REDUCTION 

2 .1. Swift Data 

A full list of the Swift observations of IC 3599 is given 
in Tabled] The Swift X-ray telescope (XRT) was operat¬ 
ing in photon counting mode and the data were reduced 
by the task xrtpipeline version 0.12.6., which is included 
in the HEASOFT package 6.12. Source counts were se¬ 
lected in a circle with a radius of 70.7 for the 2010 Febru¬ 
ary and May observations, and 47.1 for all others. The 
larger source extraction radius during the 2010 
observations was chosen because the source was 
brighter. The background counts were collected in a 
nearby circular region with a radius of 247.5 . Only the 
2010 February and May data had enough counts to al¬ 
low for a spectral analysis of each individual observation. 
The data of 2013 and 2014 were merged to allow for a 
spectral analysis. For all spectra we used the most re¬ 
cent response file swxpc0tol2s6-20130101v014-rmf and 
XSPEC version 12.7.1 for analysis. Due to the small 
number of counts used in the spectra, the counts were 
not binned. 

The UV-optical telescope (UVOT) data of each seg¬ 
ment were coadded in each filter with the UVOT task 
uvotimsum. Except for 2013 November, all observations 
were performed in all 6 UVOT filters. Source counts in 
all 6 UVOT filters were selected in a circle with a ra¬ 
dius of 3 and background counts in an annulus around 
the source extraction region, with radii of i?in = 5 
and i?out = 20 , in order to subtract the contribu¬ 
tion of star light from the host galaxy. UVOT mag¬ 
nitudes and fluxes were measured with the task uvot- 
source and aperture corrected by setting the parameter 
apercorr=CURVEOFGROWTH based on the most re- 
cent U VOT calibration as described in iBreeveld et al.l 
(Imoll . The UVOT data were corrected for Galactic 
reddening (En = 0.015; ISchleeel et al.l ITO^ . using 
equation (2) in IRoming et al.l (|2009ll . 

2 .2. XMM-Newton Slew Survey 

The field of IG 3599 was covered during t he XMM- 
Newton slew 9137800003 (|Saxton et al. Il2008ll . on 2007- 
June-19. No photons from the source were detected and 
the EPIC pn, 2-sigma, upper limits of 0.8 and 0.7 counts 
s-1, imply flux upper limits of 1.2x10“^® and 9.6x10“^® 
W m^ in the 0.2-12 keV and 0.2-2 keV bands respectively, 
assuming a spectrum with slope ax=2.5 and Galactic 
absorption of Vh = 1.16 x 10^*^ cm“^. 

3. RESULTS 

3.1. X-ray Spectral Analysis 

The X-ray spectrum during the outburst in 2010 
February can be fitted with a simple absorbed power 
law model with a X-ray spectral index ax=2.63±0.30 
and the absorption column density fixed to the Galactic 
value. The flux in the 0.2-2.0 keV band assuming this 
model is (1.41±0.11)xl0“^'^ W m“^. This is about half 
the flux that was detected during the RASS (IGruoe et al.l 
The 0.3-10 keV flux was 9.26;^ j^ x 
W m“^. The spectral slope is similar to that measured 
during the RASS. Alternatively, the X-ray spectrum can 


also be fitted by an absorbed blackbody plus power law 
model (with ax fixed to 1.0) which results in a black- 
body temperature equivalent of kT = 91±16 eV. This is 
the same value as measured during th e original outburst 
during the RASS (IGruoe et al.lll99^ . 

The Seyfert 1.9 classification of IC 3599 does suggest 
that we should expect to see excess absorption above the 
Galactic value. Although there is some hint of excess 
absorption when the spectrum is fitted by a single power 
law model, there is no evidence for excess absorption 
when the data are fitted by a blackbody model. How¬ 
ever, the quality of the relatively short-exposed spectrum 
during the high state in February 2010 does not allow fi¬ 
nal conclusions on the intrinsic column density. 

In order to obtain an X-ray low state spectrum, we 
combined the data from the Swift observations of 2013 
and 2014 which resulted in a total exposure time of 16.9 
ks with a total of 50 source counts. The s pectral analy - 
sis was performed applying Cash statistics (iCash Ill9^ . 
This spectrum can be modeled with a harder spectrum 
with an X-ray spectral index of ax=l-40±0.70 with Ah 
fixed to the Galactic value. There is no indication 
of excess absorption beyond that. The 0.3-10 keV 
flux of this combined spectrum is (9.67±1.62) x 10“^^ W 
m“^ - a factor 100 lower than during the 2010 February 
outburst. Note that when applying the absorbed black¬ 
body plus power law model mentioned above during the 
high state, the black body component has disappeared 
in the low state spectrum. 

3.2. Long-term Light Curve 

The fluxes in X-rays and in the UV/optical obtained 
by Swift are listed in Tabled] Th e ROSAT ol j serva- 
tions have already been p ublish ed in IGruoe et al.l (jl995l 
1200111 and iBrandt et al.l (1199511 . and the C h andra mea- 
surements were discussed in I Vaughan et all (|2004ll . The 
long-term light curve including our latest Swift obser- 
vations is shown in the left panel of Figure d] (see also 
iKomossa et al.l 120151: iCamoana et al.l 1201511 . together 
with the XMM-Newton slew-survey upper limit. 

3.3. Catalina Sky Survey Optical Light Curve 

While we do not have X-ray observations between the 
Chandra and XMM-Newton slew survey observations in 
2002 and 2007, respectively, IC 3599 wa s monitored 
in th e optical by the Catalina sky survey (iDrake et al.l 
1200911 . The Catalina light curve is shown together with 
the Swift XRT and UVOT W2 light curves in the right 
panel of Figured] According to the Catalina light curve, 
IC 3599 was already optically bright in 2008, and the 
outburst must have happened around MJD 54570 which 
translates to 2008 April 14. Then, either a second optical 
outburst occurred in 2010, or else the outburst lasted for 
roughly two years with the peak at around MJD 55260 
(2010 March 05) which was about a week after the Swift 
observation on 2010 February 25. These observations 
provide potentially tight constraints on outburst scenar¬ 
ios. 

3.4. Spectral Energy Distributions 

The spectral energy distributions during the February 
2010 high state and the 2013/2014 low state are displayed 
in Figured] The drop seen in X-rays is significantly larger 















































A new outburst of IC 3599 


3 


than in the optical/UV (note, however, that the optical 
emission may still be affected by host galaxy contribu¬ 
tion. If so, these values are upper limits for the core 
emission). Under these constraints, and assuming a black 
hole mass of 1 O^M 0 (see below), the Eddington ratios 
during these two epochs were L/LEdd=0.04 and 0.003, 
respectively. 

3.5. The Optical SDSS Spectrum 

A new optical spectrum of IC 3599 was t aken in the 
cours e of the Sloan Digital Sky Survey ISDSS. lYork et al.1 
1200011 in 2005-December-25 (MJD 53729). We have used 
this spectrum for a SMBH mass estimate, and in order 
to check, whether e mission lines have flared again, as 
previously observed (jBrandt et al.lll99'^ . We find that 
narrow emission lines including [FeVII] transitions are at 
their low emission levels, as seen in previous low-state 
optical spectra from 1995 (|GruDe et al.lll99^ and 1997 
(|Komossa fc Bade! 1199911 , and emission-line ratios are 
consistent with little or no variability since then. An 
analysis of the Ha emission line complex shows, that the 
spectrum is still consistent with the presence of a faint 
broad emission line component, with FWHM(Ha)^ 1020 
km s“^, similar to the value measured before (Brandt 
et al. 1995, Grupe et al. 1995, Komossa & Bade 
1999). Because of the Seyfert 1.9 classification, we do 
not use this line width to estimate the SMBH mass of 
IC 3599, since application of the well-establish ed width 

luminosity relation (e.g. iPeterson et al.l[2004ll requires 
an unobscured view on the bulk of the broad-line re¬ 
gion (BLR). Instead, we use the width of [OHI A500 7 
as a proxy for velocity dispersion (e.g. iNelsonI 12000( 1. 
and then use the scaling relati on with SMBH mass as 
given bv lTremaine et all (1200211 . First, we have fit the 
[OHI] emission with a Gauss profile, and have care¬ 
fully checked that [OHI] is well represented by a sin¬ 
gle component, and does not show a strong asymmetry, 
or an extra blue wing. We then obtain a line width of 
FWHM([OHI]) = 280 km/s (corrected for instrumental 
resolution), well consistent with the value of 260 km/s 
reported bviKomossa &: Bade I (11999(1 . Applying the re¬ 
lation of iTremaine et al.l (|2n02[ l then gives Mbh = 1-2 
xlO^ Mq. 

4. DISCUSSION 

We have detected, and followed up, a second outburst 
in the Seyfert galaxy IG 3599 with Swift. This has impor¬ 
tant implications for outburst scenarios of this unusual 
source, because some of those initially considered, did 
not predict repeat flaring. Since the value is important 
for the further discussion, we first comment on several es¬ 
timates of the SMBH mass of IC 3599, and then proceed 
with discussing outburst scenarios. 

4.1. Black Hole Mass Determination 

We have followed several different approaches of es¬ 
timating the SMBH mass of IC3599. First, applying 
scaling relations between velocity dispersion crjoiii] and 
SMBH mass, we obtain Mbh = 1-2 xIO^Mq (Sect. 
3.5). Second, from the relat ionship of black hole m ass 
to bulge K-band luminosity ([Marconi fc Hunt! 1200311 we 
find Mbh ~ 2.2 x 1 O®M 0 , with a systematic error of 0.3 


dex, from the 2MASX magnitude of fc = 11.97, after sub¬ 
tracting 25% contributed by the central AGN-dominated 
point source and co rrecting for the bulge ratio of 0.42 
(iSimard et al. l[ 2 ?)TTI l. Third, an order of magnitude es¬ 
timate of the SMBH mass can also be derived from the 
temperature of the accretion disk measured from the 
high-state spec trum in 2010, kT = 91 eV. Following 
iPetersorJ (119971 equation 3.20), weobtain~ 2 . 5 xlO®M 0 . 

All three estimates are in the range (2-12) x 10® M 0 . 
We note in p assing, that this is mu ch higher than the 
value used by iGampana et al.l (l2(315ll. 3 x 1 0 ®Mf 7 >, who 
used the width of narrow H/3 of iSani et al.l (1201(111 and 
applied BLR scaling relations. 

4.2. Outburst Scenarios 
4.2.1. Repeat Tidal Disruption Events 

At first glance, the detection of repeated flaring is un¬ 
expected, if the TDE interpretation of IC 3599 was cor¬ 
rect, because classical TDEs are rare events , at a rate 
of about 10~^-10~® p er year and per galaxy ([Alexandeil 
l2012t iKomoss^ 120121 for reviews of theory and obser¬ 
vations, respectively). However, recent work has shown 
that recurrent outbursts or high states may also occur 
in TDE scenarios, for instance when one of the follow¬ 
ing conditions is met: (1) If IC3599 hosts a binary 
black hole or recoiling black hole, tida l disruption rates 
are temporarily strongly boosted (e.g. iChen et al.l 1^0091 
iKomossa fc MerritFI 12(1081 iStone fc LoebI 1201111 . and a 
new disruption event can occur within decades. (2) If 
a TDE in IC 3599 happened in a binary SMBH, its light 
curve would show characteristic recurrent dips, since the 
presence of the secondary tempo rarily interrupts the ac¬ 
cretion stream on the primary (|Liu et al.l 120141 1. How¬ 
ever, in both cases it is unlikely that they produce a peak 
luminosity very similar to the e arlier one in 1990 . (3) A 
third possibility was studied bv iCampana et afl (I2015I1 . 
who suggested that the flaring is due to repeat tidal strip¬ 
ping of an orbiting star. However, this scenario, which 
predicts a rapid rise to maximum (their Fig. 2) may be 
difficult to reconcile with the optical high-state already 
seen in 2008. 

4.2.2. Variability Related to a Supermassive Binary Black 

Hole 

If the outbursts of IC 3599 continue repeating, the 
behavior is reminiscent of the blazar OJ 287 (e.g. 
iSillanpaa et al.llT9^ iValtonen et al.ll2008ll . and a binary 
SMBH might be responsible; i.e., a secondary BH im¬ 
pacting the accretion disk around the primary while or¬ 
biting. The two optical high-states in 2008 and 2010 
are then similar to the double-peaked optical high-states 
of OJ 287, which have been interpreted as the two disk 
impacts of the secondary BH during its orbit. Another 
mechanism which produces repeat outbursts is episodic 
stream-feeding of one of the black holes in a binary 
SMBH system (|Tanakall2013fl which recurs every orbit. 
Ongoing monitoring with Swift is required in order to 
test these scenarios further. 

4.2.3. A Disc Instability 

Given the evidence that IC 3599 hosts a long-lived 
AGN, extreme processes in its accretion disk are a 
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likely cause of the repeat outbursts. These are sug¬ 
gestive of a recurrent regular process. Similar large 
flares have been seen in a small number of accreting 
Galactic sources (e.g. GRS 1915-1-105, GRO J1744- 
28, a r id IGR 17091-3624: iBelloni et al.l[T997at iGannizz^ 
I199(H iPahari et al. 1 1201, H respectively). In these cases 
the fast rise and decay in flux has been attributed to 
an i nstability in a radiation-pr essure do minated inner 
disk (iLightman fc Eardlevll974ll . Recently. iSaxton et al.l 
(1201,511 have applied the instability to a single flare seen in 
the nucleus of the galaxy NGG 3599. In this mechanism, 
the inner disc is initially empty or filled with tenuous gas 
out to a truncation radius, Rtrunc- During the quiescent 
phase the inner disc is filled by diffusion from the outer 
disc without significantly increasing the X-ray flux. At 
some point, the internal radiation pressure of the disc 
becomes greater than the gas pressure causing a heating 
wave, that increases the local viscosity, scale height and 
accretion rate, to propagate from the inner radius back 
to the truncation radius. In this phase the luminosity 
increases quickly. Later, as the inner disc is accreted 
more rapidly than it is replenished, the flare decays back 
to the quiescent level, a new hole is left in the central 
region and the cycle begins again. 

The repeat tim e for flares is given by the filling time 
of the inner disc (|Saxton et al.ll20I5ll as 


in units of months where Mg is the black hole mass 
in units of 10 ® Mq and Rg = GM/c^, the gravitational 
radius. Assuming that the infill rate is Eddington limited 
{Medd = l)j the viscosity a = O.I and the inner radius 
Rq = 3Rg then for a flare repeat time of 19.5 years and 
Mbh = 10® — IO^Mq we find i?tr«nc=5-45 Rg. This 
is comparable to the value of 4-22 Rg ( depending on 
spin) w hich was inferred in GRS I9I5-I-I05 (|Belloni et al.l 
II997airgi . 

A lower limit for the rise time i s given by the time 
taken to fully hea t the inner disc (IBelloni et al.l[T997at 
ISaxton et al.ll20I^ : 


Trise ^ 1-5 X 10^ ^ Rt^nc ^ ^2) 

For the flare in IC 3599, r^se ^7 — 8 days, consistent 
with the Trise < 2.4 — 17 months inferred from the 
Catalina curve. 0 This scenario is supported by the 
X-ray spectra which can be explained by a weak power- 
law component augmented by a highly variable thermal 
component. 


Tfill 


0.33a-®/i°Mg/'^M”y^° 


R. 


trunc J 
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Rq 
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TABLE 1 

Swift XRT AND UVOT fluxes of IC 3599 


Tstart ^ 

MJD 

Fo.3-10keV^ 

Fv® 

Fb® 

Fu® 

Fwi® 

Fm2® 

Fw2® 

2010-02-25 07:46 

55252.365 

{9.26t°oit) X 10"^" 

8.19±0.40 

6.72±0.33 

5.64±0.34 

4.73±0.40 

4.63±0.28 

4.10±0.35 

2010-05-17 00:58 

55333.049 

(4.77±0.82)xl0-i® 

7.39±0.46 

6.11±0.33 

5.01±0.37 

4.17±0.43 

4.39±0.31 

3.95±0.35 

2013-10-30 00:34 

56595.160 

(1.051° ®?) X 10-16 

4.67±0.27 

3.85±0.20 

2.04±0.15 

1.54±0.14 

1.26±0.10 

1.31±0.11 

2013-11-06 02:24 

56602.208 

(l.Oll? ") X 10-16 





1.26±0.07 


2014-03-26 07:50 

56742.583 

(7.12l?-64) X 10-1^ 

4.94±0.39 

4.11±0.27 

2.24±0.20 

1.56±0.18 

1.44±0.09 

1.36±0.14 

2014-08-08 07:07 

56877.306 

< 2.9 X 10-16 

4.91±0.37 

4.00±0.25 

2.21±0.19 

1.73±0.18 

1.47±0.14 

1.33±0.13 

2014-11-15 04:12 

56976.188 

(l.5ll°;f6) X 10-16 

5.03±0.42 

4.02±0.28 

2.03±0.12 

1.38±0.16 

1.14±0.13 

1.20±0.13 

2014-11-23 08:57 

56984.375 

< 6.9 X 10-16 



2.10±0.15 




2015-03-29 21:37 

57110.969 

(4.7l®;») X 10-17 

4.84±0.28 

3.75±0.21 

2.02±0.15 

1.32±0.13 

1.20±0.10 

1.08±0.10 


^ Start times are given in UT 

^ Observed 0.3-10 keV flux in units of W m~^. Fluxes from the observations of 2013 and 2015 were based on an energy conversio n 
factor derived from the 2010-05-17 data and the uncerta inties were determined using the method described in IKraft et al. I Ill991i) . 
Upper limits are 3a as described in IKraft et al. 1 11199lib 

^ The reddening corrected fluxes in the UVOT filters are given in units of 10“^^ W m“^. 



Fig. 1.— Long-term 0.2-2.0 keV light curve of IC 3599, starting with the RASS observation. 

Vaughan, S., Edelson, R., & Warwick, R.S., 2004, MNRAS, 349, Wright, E.L., 2006, PASP, 118, 1711 

LI York, D.G., Adelman, J., Anderson, J.E., et al., 2000, AJ, 120, 

Valtonen, M., Lehto, H.J., Nilsson, K., et al. 2008, Nature^ 452, 1579 

851 
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Fig. 2.— Swift XRT and W2 light curves of IC 3599 and in the lower panel the optical light curve obtained by the Catalina Sky Survey. 
XRT flux is in units of 10“^® W and W2 flux is in units of 10“^^ W m~^. 



Fig. 3.— Development of the spectral energy distribution between the 2010 February high state and the 2013/2014 low state. 









